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DEFINITION OF SYMBOLS 


Symbol 



c 


m 

n 

m 

n 


Definition 
matrix coefficients 

distance from ring neutral axis to corrugation centerline, 
positive when rings are on outside 

number of longitudinal half waves 

number of circumferential waves 

= m 7r/L 
= n/R 


q 


buckling load per unit length of circumference, measured 
at corrugation centerline 


t 


corrugation thickness 


t , t , 
x y 


t 

s 


u, V, w 


X, y, z 


equivalent shell thicknesses in the x and y directions and 
for shear, respectively per unit length of circumference 

displacements of cylinder at corrugation centerline 
inx, y, and z directions, respectively 

coordinates in axial, tangential and radial directions, 
respectively 


A 

r 


area of ring 


D 

x 


El 

x 


D , D 
y xy 


circumferential bending stiffness and torsional stiffness of 
the corrugated cylinder, respectively 


D 

yr 

D 

zr 

E, E 

r 

E 


= E I 
r yr 


/L 


r 


= E I /L 


r zr r 

moduli of elasticity of corrugation and ring, respectively 


= E t 


x 


= E A /L 
r r r 
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DEFINITION OF SYMBOLS (Cont'd) 


Symbol 
G , 


I 

X 


I 

yr 

I 

zr 

J 

r 

K 

r 

L 


L 

r 


M , M 
x y 


M 

xy 

M , M 
yr zr 


M 

yxr 

N , N 
x y 


N , N 
xy yx 

N 

yr 

N 

yxr 


Q x’ Q y 


Q 


yr 


Definition 

shear moduli of corrugation and ring, respectively 
= Gt 

s 

moment of inertia of corrugation per unit length of circum- 
ference 

moment of inertia of ring about its centroid in plane of curvature 

moment of inertia of ring normal to plane of curvature 

l/G times torsional stiffness of ring 
r 

= G J /L 
r r r 

length of cylinder 
ring spacing 

stress couples acting on skin element in x and y directions, 
respectively 

torsional stress couple acting on skin element 

stress couples acting on ring element in y and z directions, 
respectively 

torsional stress couple acting on ring element 

stress resultants acting on skin element in x and y directions, 
respectively 

shear stress resultants acting on skin element 

stress resultant acting on ring element in y direction 

shear stress resultant acting on ring element in x 
direction 

radial shear stress resultants acting on skin element 
radial shear stress resultant acting on ring element 
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DEFINITION OF SYMBOLS (Concluded) 


Symbol Definition 

R radius to centerline of corrugated skin 


R 

r 

T , T 
x y 

u, v, w 
x, Y, Z 

a 

P 


= R + c 

moments acting at interface of skin and ring, per unit area 

amplitudes of u, v, and w, respectively 

forces acting at interface of skin and ring, per unit area 

number of rings per longitudinal half wave length 

bending factor, ratio of maximum buckling stress in pure 
bending to buckling stress in axial compression for a given 
cylinder 

ratio of maximum bending stress to maximum total stress 
for a cylinder under combined loading 

shear strain of corrugation 

strains of corrugation in x and y directions, respectively 
circumferential strain of ring 

= y/R 

curvature changes of corrugation in x and y directions, 
respectively 

specific twist of skin element 

curvature changes in plane and normal to plane of ring, 
respectively 

specific twist of ring 

= q m 2 

Poisson’s ratios associated with bending in the x and y 
directions, respectively 

Poisson’s ratios associated with extension in the x and y 
directions, respectively 

longitudinal stress in corrugation 
vii 
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THE GENERAL INSTABILITY OF RING-STIFFENED 
CORRUGATED CYLINDERS UNDER AXIAL COMPRESS ION 

SUMMARY 


A method is presented to determine the general instability load of a ring- 
stiffened corrugated cylinder under axial compression. This method is developed 
using linear small deflection theory. The stiffness properties of the rings are 
uniformly distributed along the cylinder and the eccentricity of the rings with re- 
spect to the corrugation centerline is taken into account. 

Analytical and experimental results are compared. In this comparison 
good agreement is obtained for cylinders loaded in pure compression. For the 
cylinders subjected to bending or a combination of bending and compression the 
analytical calculations are conservative. 

A computer program for the application of this method has been developed. 
The program and an explanation of its notation are included in this report. 


INTRODUCTION 


Axially loaded ring-stiffened corrugated cylinders are susceptible to two 
types of instability; panel instability and general instability. Panel instability is 
defined as buckling of the cylinder between rings. In general instability, the rings 
as well as the cylinder undergo buckling displacements. Panel instability is dis- 
cussed in many textbooks and manuals; however, the same is not true for general 
instability. 

In comparing the stiffness characteristics of a corrugated cylinder to those 
of a monocoque cylinder of the same weight it must be noted that the corrugated 
shell has a much greater longitudinal bending stiffness. This advantage is offset 
to some degree by the low circumferential extensional stiffness of the corrugated 
cylinder. Because of this lack of circumferential extensional stiffness each 
corrugation will act as an individual column unless the shell is reinforced by rings. 
Low extensional stiffness may be desirable, however, in areas of high thermal 
gradients. 

In the past, ring- stiffened cylinders have often been designed to prevent 
general instability by employing semi-empirical methods to size the rings. One 
of the most commonly used of such methods is given by Shanley [ 1] . Until re- 
cently it was believed that the ring moment of inertia obtained from his formula 



was sufficient to safeguard against general instability. Tests performed under 
the direction of the Marshall Space Flight Center as part of the Saturn V develop- 
ment program have shown, however, that Shanley’s criterion may be very uncon- 
servative for the design of ring-stiffened corrugated cylinders. The results of 
these tests are given in Tables I and II. 

Since Shanley's work was published, considerable research has been 
performed on the stability of stiffened cylinders. A major contribution was made 
by Van der Neut [ 2] , when he considered the effect of ring and stringer eccentric- 
ities. Additional contributions were made by Hedgepeth and Hall [3] , Card [4], 
anrl Baruch and Singer [5] . Often, however, research in this field has produced 
a method which is either too academic or too complex to be used by the stress 
analyst. 

The purpose of this report is to present a reliable and relatively simple 
means of predicting the general instability load of corrugated cylinders under axial 
compression. This method considers the eccentricity of the rings with respect 
to the corrugation centerline. It also incorporates all the stiffnesses attributable 
to the rings and the shear and longitudinal stiffness properties of the corrugation. 
The circumferential extensional and bending stiffnesses and the torsional stiffness 
of the corrugation are small (generally less than 1 percent of the longitudinal 
stiffnesses) and are therefore neglected. 


GENERAL THEORY 
Assumptions 

The method of analysis presented in this paper is based on the following 
assumptions: 

1. Linear small deflection theory applies. This is justified because 
the high longitudinal bending stiffness of a corrugated cylinder makes it less 
susceptible to initial imperfections and other monocoque effects. 

2. The longitudinal wavelength of the buckled skin is sufficiently large 
to permit a "smeared" ring approach; i. e. , all the ring stiffness parameters 
may be uniformly distributed along the cylinder. 

3. The corrugated cylinder can be treated as an equivalent orthotropic 
cylinder, the radius of which is equal to the mean radius of the corrugation. 
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4. The circumferential extensional and bending stiffnesses of the cor- 
rugation as well as its torsional stiffness can be neglected since they are small 
as compared to the longitudinal stiffness properties of the shell. 

5. Buckling displacements are sinusoidal. 

6. Prebuckling deformations are neglected. 

7. Plasticity effects and local failures are not considered. 

Displacements and Boundary Conditions 

The cylinder is in equilibrium under the applied load just prior to buckling 
and deformations due to buckling are measured from this position. In accordance 
with assumption 5, displacements may be written in the form 


TT H17TX „ 

u = U cos — — cos n0 


TT . m7rx . 

v = Vsm- — sm no (1) 

Li 

„ T . m7TX 

w= W sin — - — cos n6> 


This corresponds to the following simply supported boundary conditions at 
x = 0, L 

w = 0 N = 0 

x 

v = 0 M = 0 

x 

Thus, at the ends of the cylinder motion radially and tangentially is prevented, 
while longitudinal motion is allowed; i.e. , u ^ 0. These boundary conditions are 
appropriate for cylinders bounded by deep supporting rings, which are rigid in 
their own planes but may readily bend or warp out of their planes. 


Elastic Relations 


For an orthotropic shell the relations for the stress resultants and couples 
in terms of the strains and curvatures are 


3 



Et 


N = 
x 


X 


x y 
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xy 
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D 
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i “M M 
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X y y 
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D 


1 jLt 

x y 
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(1 K ) 
X X 


(2) 


M = D k 
xy xy xy 


In view of assumption (4) and with the use of the reciprocal theorem, 
equations (2) lead directly to the following expressions for the corrugated cylinder 


N = E € 

X X 



N 

xy 


= G ( 


R 


3u 

90 



M = D 
x x 


K 

X 


= D ts. 

x 9x 2 


(3) 


N = M =M = 0 
y y xy 

The rings are considered eccentric with respect to the centerline of the 
skin and displacements of the ring at its centroid may be found by means of the 
transformations 


u = u 
r 


c 


9w 

9x 


v 

r 


_r _ jc 9w 
R v “ R de 


( 4 ) 


w 

r 


= w 
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where c is the distance from the ring centroid to the skin centerline. For a ring 
loaded normal to, as well as in the plane of, initial curvature, the strain, curva- 
tures and specific twist at the central axis, may be written 


= x 

e yr R 

r 



w 

r 


K 

yr 




K 

zr 




( 5 ) 


K 

yxr 




a 2 


w 


8x dO do 


Assuming the shear center to coincide with the ring centroid, the stress 
resultants and stress couples acting on an element of the ring are related to the 
strains and curvatures by the equations 


N = E c 
yr r yr 


M = 

D 

K 

yr 

yr 

yr 

M 

D 

K 

zr 

zr 

zr 

M 

- K 

K 

yxr 

r 

yxr 


Substituting equations (4) and (5) in the above expressions gives 


N ~ E 
yr r 


1_ 

R 


dv 1 
dO + R 

r 



a 2 w 
do 2 


M 

yr 


D 

yr 

R R 

r 


f j3v , 
do do 2 


M 

zr 



a 2 u _ _a_ 

dO 2 dx 



( 6 ) 
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(6 concluded) 
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yxr 
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r 




cfiw A 
dx 36 J 


Equilibrium Equations 


Differential equations are obtained by considering separately the equi- 
librium of an element of the skin and a corresponding "smeared" ring element. 

Six conditions of equilibrium must be satisfied by the forces and moments acting 
on each of the elements. An element of the skin is shown in Figure 1 a, b, c. 

The first of these figures shows the stress resultants and the forces X, Y, and Z 
acting at the interface of the skin and the ring, and the second figure gives the 
stress couples and the moments T x and Ty transferred into the skin element from 
the ring. Components of the buckling force resulting from the deformation of the 
element are shown in Figure lc. The component in the x direction is multiplied 
by the quantity (1 + e x ) to include the straining of the middle surface as suggested 
by Flugge (7) . The components in the y and z directions are due to the rotations 
9v/ dx and 9w/9x of the element respectively. 

The six conditions of equilibrium for the skin lead to the following equations 


9N , 9N 
x 1 yx 

9x R 90 
9N 

sr 

9x 

9Q 
x_ 

ax 


8M 
x 

dx 


+ Q 


N - N 
xy yx 


q|| + X=0 

q§ +Z =° 


+ T = 0 

x 


T = 0 

y 

= o 


(7) 


These equations may be combined and put in the form 
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(8) 


9N , 9N 

x i yx . 9ni 

= — — J — i- q — 5- 

9x R 99 4 9x 2 
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R 

9T 

x 
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9N *2. 

IS afv 

9x ^ 9x 2 


J. 
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9T 

S 
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9 2 M 
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+ q 


d^w 

9x 2 


The six equations of equilibrium for an element of the ring (Fig. 2) may 
be obtained in a similar manner. They are 
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The last three of the above equations are used to eliminate Qy r and Ny xr - 
This yields the following three equations 
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( 10 ) 
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The forces X, Y, and Z and the moments T x and Ty acting at the inter- 
face of the skin and ring may now be eliminated by subtracting equation (8) from 
equation (10) . Hence, 


9N , 9N . S 2 M , 9M 

x + _1 yx _ 1 zr 1 yxr 

9x R 90 RE 90 2 R R 80 
r r 


9\i 


3N R 9N , 9M 
yx yr 1 yr 

9x R 2 90 ~ R 2 de 

a 2 M n a 2 n . a 2 m 
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aV 
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, a 2 m 

1 yxr 

r dx de 
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J. 
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8M 

zr 

dx 


c_ 

R R 


a 3 m 

zr 


8x 90 2 q 


a 2 w 

9x 2 


= 0 


( 11 ) 


r 

These are the equilibrium equations in terms of the stress resultants. 
Substituting the expressions (3) and (6) for the stress resultants into equations 
(11) gives the equilibrium equations in terms of the displacements. 
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R R 3 
r 



+ R 


a 3 w A 
9x de 2 J 


- q 


a^u 

9x 2 


= 0 


— d 2 u , aVA — f ^r aV l 8w c a 3 w ^ 
R 9x 90 dx 2 ] r VR 3 00 2 R 2 de ~ R 3 dd 3 J 


+ 


D 

yr 

R 3 R 

r 




= 0 


-D 

x 


a 4 w , E r 
8x 4 R 2 


dv c 9 3 v R 2c 8 2 w 
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c 2 a 4 w ^ 

R R 90 4 J 
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( 12 ) 


^yr f a 3 v _ a 4 w A _ ^r f a 3 u 8 4 w ^ ^zr f a 3 u 

R 3 R r ^90 3 " 80 4 J ~ R 3 \8x80 2 + H 9x 2 90 2 J “ R R 2 ^8x 80 2 + 


c a 5 u a^ a 4 w c 2 a 6 w A a^ 

R dx 80 4 " r dx 2 - ° dx 2 de 2 ~ R 8x 2 80 4 / “ q 9x 2 ~ 
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After introducing the expression for the displacements given by equations 
(1) into the differential equations (12) one obtains the following three linear equa- 
tions . 



+ q m 2 W = 0 


Determination of Buckling Load 

The set of homogeneous equations (13) has nontrivial solutions only when 
the determinant of its matrix is zero. Equations (13) may be written in matrix 
form 


a ll + ^ 

a 12 

a 13 

a 21 

a 22 + ^ 

a 23 

a 31 

a 32 

a 33 + X 


(14) 
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The determinant of the characteristic matrix of A is a polynominal of the 
third degree in A. The problem of determining the buckling load (q = A/m 2 ) for 
known values of m and n has therefore been reduced to that of finding the roots 
(eigenvalues) of the characteristic equation 

a + bA + cA 2 + A 2 = 0 (15) 

where a, b, and c are known functions of the coefficients a^j. Obviously only 
real and positive roots of equation (15) are of interest, the lowest of which de- 
termines the buckling load for the mode shape under consideration. The critical 
buckling load of the cylinder may be found by calculating q for a range of values 
of m and n, and plotting a family of curves as shown in Figure 3. The critical 
buckling load will then be the minimum value of q corresponding to integer values 
of m and n. If the computer program (see Appendix) is used, the minimum 
buckling load will be indicated for the specified range of m and n, and no plotting 
is required. 

Since W is usually in the order of n times larger than V or U an approxi- 
mate solution for q may be obtained by dropping the terms containing q in the 
first two of equations (13). This yields 


q m,n (a 2 2 - an a 22 ) m 2 

where |A| is the determinant of the symmetrical matrix 


ail 

a 12 

a i3 

a 12 

a 22 

a 23 

ai3 

a 23 

a 33 


(16) 


If an iterative procedure is used to find the roots of equation (15) , the above 
approximation is found to be very useful as a starting point. 

COMPARISON WITH TEST RESULTS 
Axial Load 


As a part of the Saturn V development program eleven ring- stiffened 
corrugated cylinders were tested to failure under axial load. Of these, five 


failed in general instability. The remaining cylinders failed either by panel in- 
stability or local crippling before the general instability load was reached. A 
photograph of a typical general instability failure is shown in Figure 4. Figure 5 
shows one of the rings of this cylinder after the failure. Table I gives the geom- 
etries of these five cylinders along with their actual and calculated failure loads. 
The close agreement between the experimental data and calculations supports the 
validity of the method of analysis given in this report. 

TABLE I. RING-STIFFENED CORRUGATED CYLINDERS - AXIAL LOADING 


CYLINDER 

No. i 

No. 2 

No. 3 

No. 4 

No. 5 

Aluminum Alloy 

7075-T6 

7075-T6 

7075-T6 

7075-T6 

7075-T6 

Cylinder Length (in. ) 

33.0 

33. 0 

33.0 

69.6 

69.6 

(cm) 

83.8 

83.8 

83.8 

176.8 

176. 8 

Radius (in. ) 

24.7 

24.7 

24. 7 

49.4 

49.4 

(cm) 

62.7 

62. 7 

62.7 

125. 0 

125. 0 

Corrugation Pitch (in. ) 

1.43 

1.43 

1.43 

2.85 

2.85 

(cm) 

3. 63 

3.63 

3. 63 

7.24 

7.24 

Corrugation Thickness (in.) 

0. 020 

0. 019 

0. 025 

0.041 

0.041 

(cm) 

0. 05 

0. 048 

a 06 

0.10 

0. 10 

Corrugation Depth (in. ) 

0.44 

0. 44 

0. 44 

0. 87 

0.87 

(cm) 

1.12 

1. 12 

1. 12 

2. 21 

2.21' 

Type of Ring 

L * 

I 

X 

X 

X 

Ring Spacing (in. ) 

6. 38 

6.38 

6.38 

12.4 

12.4 

(cm) 

16.21 1 

16.21 

16. 21 

31. 5 

31.5 

Ring Moment of Inertia (in. 4 ) 

0. 0050 , 

0.0104 , 

0. 0104 

0. 286 

0.286 

(cm 4 ) 

0. 208 

0. 433 

0. 433 

1 

11. 9 

11.9 

Ring Area (in. 2 ) 

0. 040 

0. 121 

0. 121 

0. 180 

0. 180 

(cm 2 ) 

0. 258 

0.78 | 

0.78 

1.16 

1. 16 

Ring Eccentricity (in. ) 

-0. 73 ! 

-0. 53 

-0. 53 

-1. 99 

-1. 99 

(cm) 

-1.85 

-1.35 

-1. 35 

-5. 05 

-5. 05 

Actual Failure Load (Kips) 

131. 

174. 

224. 

659. 

648. 

(N) 

5. 83 x 10 5 

7 . 74 x 10 5 ! 

9. 96 x 10 5 

2.93 x 10 G 

2. 88 x 10 6 

Calculated Failure Load (Kips) 

118. 

198. 

233. 

654. 

654. 

(N) 

5.25 x 10 5 

8. 81 x 10 5 

1 . 04 x 10 6 

2.91 x 10 5 

2.91 X 10 6 

Percent Error (%) 

11.0 

-12.1 

-3. 9 

0.8 

-0.9 


* For this ring an effective 
moment of inertia and ef- 
fective area were calculated 
using the approach given in 
Reference 9. 
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Ill 


The instrumentation on these cylinders provided some interesting informa- 
tion. The radial deflection gages located on the rings showed that as the load 
increased most of the rings deflected into the theoretically predicted circumferen- 
tial mode shape. Figure 6 shows a typical example of this phenomenon. When the 
radial deflection at a point of maximum bending on the ring is plotted versus axial 
load the curve obtained is asymptotic to the failure load as is shown in Figure 7. 

Another interesting phenomenon was shown by the strain gages located on 
the corrugated skin. Two gages were located opposite each other at various points 
on the skin, see Figure 8, so that longitudinal bending of the corrugation could be 
observed. As the axial load on the cylinder increased, the strain versus load plot 
of the two gages in each set remained linear and almost coincident until just below 
the failure load. Then at various locations on the cylinder the strain in one gage 
increased nonlinearly while the strain in the other leveled off and then decreased as 
shown in Figure 9. This meant that the corrugation was bending appreciably, and 
failure of the cylinder was imminent. The actual cylinder failures were similar to 
monocoque failures in their swiftness. Conjecture can be made that this second 
phenomenon might be used to determine the general instability failure load of a 
ring-stiffened corrugated cylinder without actually failing the specimen. 


Axial Load and Bending 

Test data is also available for a ring-stiffened corrugated cylinder loaded 
in pure bending and for a cylinder loaded simultaneously in bending and axial com- 
pression. Table II gives the cylinder geometries and loads. Cylinders loaded in 
the above manner are usually analyzed by calculating the maximum compressive 
stress due to bending and axial load and then assuming this stress to act along the 
entire periphery of the cylinder. This, however, leads to a conservative prediction 
for the buckling stress, since a stiffened cylinder can withstand a greater maximum 
stress in bending than in pure compression. This increased load carrying capability 
of cylinders in bending can be expressed in terms of a bending factor, /3 . If the 
bending factor /3 and the allowable stress in pure compression a are known then the 
maximum allowable stress in any combination of bending and compression may be 
calculated from the equation 


cr 

max 
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TABLE II. RING-STIFFENED CORRUGATED CYLINDERS - 
AXIAL AND BENDING LOADING 


CYLINDERS 

No. 6 

No. 7 

Aluminum Alloy 

7075-T6 

7075-T6 

Cylinder Length (in. ) 

268.6 

33.0 

(cm) 

682.2 

83.8 

Radius (in. ) 

197. 6 

24.7 

(cm) 

501. 9 

62.7 

Corrugation Thickness (in. ) 

0. 145 

0.020 

(cm) 

0. 368 

0. 051 

Corrugation Pitch (in. ) 

11.40 

1.43 

(cm) 

28. 95 

3.63 

Corrugation Depth (in. ) 

3.48 

0.44 

(cm) 

8. 84 

1. 12 

Type of Ring 

1 

c* 

Ring Spacing (in. ) 

51.0 

6.38 

(cm) 

129. 54 

16.20 

Ring In- Plane Moment of Inertia (in. 4 ) 

34.4 

0. 0050 

(cm 4 ) 

0. 1431 

0.208 

Ring Out-of- Plane Moment of Inertia (in. 4 ) 

0 

0 

(cm 4 ) 

0 

0 

Ring Torsional Stiffness/ G R (in. 4 ) 

0 

0 

(cm 4 ) 

0 

0 

Ring Area (in. 2 ) 

2.48 

0.040 

(cm 2 ) 

16. 00 

a 258 

Ring Eccentricity (in. ) 

-6.24 

-0.73 

(cm) 

-15. 84 

-1. 85 

Actual Bending Load (Kip-ft. ) 

27, 900 

161.0 

(N-m) 

3. 79 x 10 7 

2. 18 x 10 5 

Actual Axial Load (Kips) 

6, 930 

0 

(N) 

30. 8 x 10 6 

0 

Actual Maximum Stress (lb. /in. 2 ) 

43, 530 

37, 960 

(N/ m 2 ) 

3.00 x 10 8 

2.62 x 10 8 

Calculated Maximum Stress without Bending Factor (lb. /in. 2 ) 

37, 560 

28, 620 

(N/ m 2 ) 

2. 59 x 10 8 

1.97 x 10 8 

Percent Error (%) 

15. 9 

32.6 

P 

1.20 

1.20 

y 

. 328 

1.0 

Calculated Maximum Stress with Bending Factor (lb. /in. 2 ) 

39, 700 

1 34, 300 

(N/ m 2 ) 

2. 74 x 10 8 

• 2.36 x 10 8 

Percent Error (% ) 

9.6 

| 10.7 

- 


L 


* For this ring an effective moment of inertia and effective area were calculated using 
the approach given in reference 9 . 


where y is the ratio of the bending stress to the maximum total ^stress . Table II 
shows the calculated maximum stress for these cylinders both without and with the 
bending factor/? . For the latter, /3 has been taken as 1. 20 which is the, ratio of the 
actual failure stress of cylinder No. 7 to the actual failure stress of cylinder No. 1. 
These cylinders are identical; one was tested in bending and the other in compression 

Additional Remarks 


The average number of rings ( a) per longitudinal half wavelength for the 
cylinders given in Tables I and II is relatively low, a varying from 1. 7 for 
cylinders No. 2 and No. 3 to 2. 6 for cylinder No. 6. Van der Neut [6 ] performed 
a study to determine what error is produced by using a "smeared" ring approach 
when a is low. He states that for stiffened cylinders the error is on the order of 
4 percent for a = 2. 0 and 6 percent for a = 1. 6, the exact error being dependent 
upon the cylinder stiffness properties. 

All the cylinders given in Tables I and II had some end fixity. Deflection 
measurements indicate that for cylinders No. 4, 5, and 6 the amount of end 
fixity was negligible. Unfortunately there were not sufficient deflection measure- 
ments on the other cylinders to determine the amount of end fixity, but it is 
believed to be small for these cylinders also. 


DISCUSSION 


The method for predicting general instability developed in this report con- 
siders the eccentricity of the ring with respect to the skin centerline. It can be 
shown that this factor has a large effect on the general instability failure load. In 
fact, moving the rings from the inside to the outside of the cylinder can sometimes 
change the general instability load 100 percent or more. As an example, the general 
instability failure load of cylinder No. 4 is 654 kips (2. 9 x 10 6 N) for inside rings 
and 1254 kips (5. 6 x 10 6 N) if the same rings are on the outside of the corrugation. 
This same type of effect is present in cylinders with inside or outside longitudinal 
stiffeners as is shown by the test data in reference 4. 

As the test results show, good agreement is obtained for the cylinders 
loaded in bending and for those loaded simultaneously in bending and axial com- 
pression, if a bending factor of i. 20 is applied to the bending portion of the load. 
Cylinders loaded in bending carry a greater maximum stress because the portion 
of the cylinder that is highly loaded is stabilized by the remainder of the cylinder. 

At the present, though, sufficient information is not available as to exactly what 
bending factor should be used for each particular cylinder, this factor being a 
function of the cylinder stiffness properties. Based on currently available infor- 
mation, it is not recommended that a bending factor be used for design. 
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The use of linear small deflection theory has been justified because the 
relatively high bending stiffness of a corrugated shell makes the cylinder less 
susceptible to initial imperfections and other monocoque effects. A corrugated 
shell for which this assumption does not hold true is not envisioned as being- 
practical, but still it would be possible to make a cylinder with such a small 
corrugation depth that it would actually approach monocoque in properties. For 
a design of this nature a reduction factor will have to be applied to the failure 
load given by this method. An acceptable technique for determining the reduction 
factor is given by Almroth [8] . 

For all the cylinders tested the ring extensional and in-plane bending stiff- 
nesses were the only ring properties affecting the failure load. It is believed 
by the authors that for most practical applications the out-of-plane bending- 
stiffness and the torsional stiffness of the ring are of secondary importance. 


In computing the stiffness properties of the ring, care should be taken by 
the analyst that the effective rather than the apparent stiffness properties are 
used. This is especially true for the out-of-plane properties such as the lateral 
bending and torsional stiffness, but may also be important for the in-plane proper- 
ties of the ring; e. g. a channel section having wide and comparatively thin flanges 
may not be fully effective in bending, as is discussed in reference 9. 

It should also be mentioned that the method given in this report only de- 
termines the general instability failure load. It does not check for panel buckling 
or local crippling failures. 


CONCLUDING REMARKS 


Linear small deflection theory has been used to develop a method to de- 
termine the general instability load of a ring-stiffened corrugated cylinder under 
axial compression. The general instability failure loads of seven corrugated cylin- 
ders have been compared with loads calculated by using this method. Agreement 
between the calculated and actual failure loads were quite close for the five 
cylinders loaded in pure compression. When a bending factor of 1. 20 was used, 
good agreement was also obtained for one cylinder loaded in pure bending 
and for one cylinder having a combination of axial and bending load. Since 
calculations must be made for many different mode shapes before the minimum 
buckling load can be determined, a computer program was developed. 
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- — Forces and moments 
assumed to be zero 



FIGURE 1. FORCES AND MOMENTS ACTING ON SKIN ELEMENT 




FIGURE 2. FORCES AND MOMENTS ACTING ON RING ELEMENT 
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FIGURE 5. RING OF CYLINDER NO. 4 



Ring Radial Deflection in. (cm) 


111 llllll 



FIGURE 6. CIRCUMFERENTIAL MODE SHAPE OF CYLINDER TESTED IN AXIAL 
COMPRESSION 
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APPENDIX 


The computer program is written in Fortran n. The program notation is 
given on page 24 and the program listing on pages 25 through 27. The sequence 
of steps for program compilation and operation is given in reference 10. 

The program input data and corresponding column locations on the input 
cards are in the following form. 



1-19 

20 - 69 

70 - 80 

Card 1 

- 

Case Title 

- 



1-16 

17-32 

33-48 

49-64 

65-80 

Card 2 

E 

ER 

G 

GR 

R 

Card 3 

QL 

TX 

TS 

QIX 

C 

Card 4 

QLR 

AR 

QIYR 

QIZR 

QJR 



1-4 

1 

00 

LO 

9-12 

13-16 

17-80 

Card 5 . 

Ml 

MM 

N1 

NN 

- 


where Ml, N1 and MM, NN are the first and last values of m and n respectively, 
and 


TX = t (D.F.) 

TS = t/D. F. 

where D. F. , the development factor, is defined as the ratio of the average skin 
area per inch of circumference to the skin thickness. 
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An example of the program output is shown on pages 28 and 29. The 
axisymmetric mode shapes (n= 0) have been included for completeness. 

In using the computer program a sufficient range of mode shapes must 
be considered so that the minimum biickling load is definitely obtained. This is 
mentioned because it is possible to have more than one apparent minimum buckling 
load as is shown in Figure 10. 

COMPUTER PROGRAM NOTATION 


TEXT NOTATION 

PROGRAM NOTATION 

a. . 

AIJ 

iJ 


c 

C 

m, n 

QM, QN 

rn 

QMBAR 

n 

QNBAR 

q 

Q3 

hi 9 

TX, TS 

A r 

AR 

D x > Dyr* ^zr 

DX, DYR, DZR 

E, E r 

E, ER 

E 

EBARX 

% 

EBARR 

G, G r 

G, GR 

G 

GBAR 

*x> ^yr» *zr 

QIX, QIYR, QIZR 

J r 

QJR 

Er 

QKR 

L, I* 

QL, QLR 

R, R r 

R, RR 
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THE GENERAL INSTABILITY OF RING STIFFENED CORRUGATED 
CYLINDERS UNDER AXIAL COMPRESSION 
1 READ 4 

READ 2 » E » ER » G» GR » R> QL » TX» TS> UIX* 

1 C, QLR » AR » QIYR, GIZR, Q JR » Ml* MM » M 1 » NN 
PUNCH’ 3 
PUNCH 4 
PUNCH 5 
PUNCH 12 
PUNCH 6 

PUNCH 7* E* ER* G» GR* R 
• PUNCH 8 

PUNCH 9* OL* TX » TS> GIX* C 
PUNCH 10 

PUNCH 9* QLR* AR* QIYR* QIZR* Q JR 
PUNCH 13 

PUNCH 14, Ml, MM* Nl» NN 

PUNCH 15 

PUNCH 16 

PUNCH 17 

PUNCH 18 

EBARX = F*TX 

GBAR = G*TS 

DX = E *Q I X 

E'BARR = ER*AR/.jLR 

DYR = ER*Q I YR/QLR 

DZR = ER*QIZR/JLR 

OKR = GR*QJR/QLK 

RR = R + C 

RRR = R/RR 

Y = 0 

20 DO 51 N = Nl* NN 
DO 51 M = Ml* MM 
QM = M 
ON = N 

QMBAR = QM*3. 14159/QL 

QNBAR = QN/R 

QR = R*C*QNBAR**2 + 1. 

QRR = -R*C*QNBAR**2 + l./RRR 

All = -EBARX *GMBAR **2 - GBAR *GNBAK**2 - ( DZ R*QNt>AR*#2 + ukR/R**2)* 
1 RRR**3*QNBAR**2 
A12 = GBAR*QMb AR*GN8AR 

A 1 3 = - (DZk*GKk + GmR) * (RRK*UNBAR ) **2*UMbAR/ RR 

A22 = - GBAR*GiiBAR**2 - ( EBARR/RRR + D YR*KRk/ R**2 ) *GNBAR * *2 

A23 = - < EBARR*GN8AR*uk/R + DYR*QNBAR**3 /RR ) 

A 3 3 = - DX*GKBAR**4 - ( EBARR*GR**2/R**2 + DYR*GNBAR**4 + RkR**2* 

1 ( DZR*QkR**2/R**2 + Q,<R *UNBAR**2 ) *GKBAR**2 ) *RRR 

DD = QM0AR**6 

CC = QMBAR** 4* ( A 1 1 + A22 + A33) 

BB = QMBAR**2* ( All *A22 + A11*A33 + A22*A3i - A1 2**2 - 
1 A 1 3 * * 2 - A2 3 **2 ) 

AA = A 1 1 * A 2 2 * A 3 3 + 2 . *A 1 2* A 1 3*A2 3 - A23**2*A11 - A13**2*A22 
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1 


- A12**2*A33 
XXX = 0. 

Q1 = 0. 

El = DD*Q 1**3 + CC*G1**2 + BB*0 1 + AA 
DQ = . 1*AA/ ( ( A12**2 - A 1 1*A2 2 ) * ( QMBAR**2 ) ) 

21 02 = Q1 + 00 

22 E 2 = DD*Q2**3 + CC*G2**2 + 66*02 + AA 
IF (E1*E2) 28 *28 *26 

26 01 = Q2 
El = E2 
GO TO 21 

28 IF (XXX) 29»29»30 

29 DO = .01*02 
XXX = 1. 

GO TO 21 

30 Q3 = Q 1 

RAD = - 3.*<DD*03>**2 - 2.*DD*CC*03 + CC**2 - 4.*DD*SB 
IF (RAD) 50. 31. 31 

31 Q4 = ( — ( CC + DD*Q3 ) + SORTF ( RAD ) ) / ( 2 . *DD ) 

Q5 = (-(CC + DD*Q3 ) - SORTF ( RAD ) ) / ( 2 . *DD ) 

IF ( Q4 ) 34. 34, 32 

32 IF ( Q3 - 04) 34. 34. 33 

33 03 = Q4 

34 IF (05) 50. 50. 35 

35 IF (03 - Q5 ) 50, 50, 36 

36 Q3 = Q5 

50 P = 6 . 28*R*Q3 
STRES = Q3/TX 

IF (Y) 55. 55. 52 
55 STREX = STRES 
MX = M 
NX = N 
GO .TO 54 

52 IF (STRES - STREX) 53. 51, 51 

53 STREX = STRES 
MX = M 

NX = N 
GO TO 51 

54 Y = 1. 

51 PUNCH 19. 03, P» STRES. M, N 
PUNCH 23 

PUNCH 24 

PUNCH 11, STREX. MX, NX 
GO TO 1 

2 FORMAT (4F16.0, F16.4/ 3 Flo. 4/ 5F16.4/ 414) 

3 FORMAT ( / / 1 9X41hTiiE GENERAL INSTABILITY UF RING STIFFEN LD/ 
H9X44HC0RRUGATED CYLINDERS UNDER AXIAL COMPRESSION//) 

4 FORMAT (19X49H ) 

5 FORMAT ( ///35X10HI NPUT DATA) 

6 FORMAT (15X1HE, 14X2HLR » 15X1HG, 14X2HGR , 15X1HR) 

7 FORMAT (4F16.0, F16.4//) 

8 FORMAT (14X2H0L, 14X2HTX » 14X2HTS » 13X3H0IX, 15X1HC) 



9 FORMAT (5F16.6//> 

10 FORMAT ( 1 3X3HQLR » 14X2HAR, 12X4HUIYR, 12X4HUIZR, 13X3H0JR) 

11 FORMAT (4XF16.4, 2(7X13)) 

12 FORMAT ( 35X10H///////////// ) 

13 FORMAT (14X2HM1* 1 4X2HMM * 14X2HN1, 14X2HNN) 

14 FORMAT (4(12X14)///) 

15 FORMAT (35X1 1H0UTPUT DATA) 

16 FORMAT ( 35X11FI / / // // /////// / ) 

17 FORMAT (15X5HAXIAL) 

18 FORMAT (7X13HLOAD PER INCH, 10X10HAXIAL LOAD, 8X12HAX I AL STRESS, 
1 9X1HM, 9X1HN/) 

19 FORMAT ( 3 ( 4X F16 • 4 ) , 2(7X13)) 

23 FORMAT ( ///3X17HMI N. AXIAL STRESS) 

24 FORMAT (2X18HIN THE ABOVE RANGE, 9X1HM » 9 X 1 H N ) 

01 FORMAT (F16.4) 

END 
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THE GENERAL INSTABILITY OF RING STIFFENED 
CORRUGATED CYLINDER'S UNDER AXIAL COMPRESSION 


CYLINDER NO. I 


INPUT DATA 
////////// 


E 

ER 

G 

GK 

R 

10500000. 

10500000. 

3900000. 

3900000. 

24.6800 

OL 

TX 

TS 

Q I X 

C 

32.950000 

♦ 0266U0 

.015100 

.000775 

-.730000 

QLR 

AR 

QI YR 

QIZR 

QJR 

6.37500C 

.040000 

.00500U 

0.000000 

0.000000 

Ml 

MM 

N 1 

NN 


1 

4 

0 

5 




OUTPUT DATA 
/////////// 


AXIAL 



LOAD PER INCH 

AXIAL LOAD 

AXIAL 8TRES6 

12333.0930 

1911820.9000 

463725.3000 

3361.1760 

52095C.0100 

126360.0000 

2028.1123 

5143.j7.9300 

7 6246- .8230 

1949.9041 

3C2216.4100 

/ 3304 • 6650 

9282.7383 

1436 ? 38 * 4000 

3469 7 :> . BbOO 

3016.3183 

46753] . 3700 

113402.9400 

1902.2971 

2948 3-7.7800 

71514.9280 

1877.6704 

291020.8800 

70589.1120 

4857.0714 

752799.4300 

182596.6600 

2262.2310 

350624.0800 

85046.2730 

1626.4831 

252069.2600 

61145.9810 

1736.1352 

269084.2800 

65268.2400 

2125.2484 

329393.0900 

79696.5560 

1445.5047 

224039.3500 

54342.2810 

1284.3015 

199054.4000 

462o2.0110 

1551.2859 

240434.4200 

58319.0130 

1219.6112 

189028.0200 

45850.0450 

889.5332 

137869.1100 

33441. 1C00 

1001.3810 

155204.4400 

37645.9020 

1386.0445 

214823.5900 

52106.9360 

1677.6855 

260025.1400 

63070.8830 

758.9283 

117626. 6000 

28531.1400 

887.8293 

137605.0100 

33377.0410 

1307.4159 

2 02636.9100 

49150.9730 


MIN. AXIAL STRESS 
IN THE ABOVE RANGE 

M 

N 


28531.1400 

2 

5 


iX 


1 

2 

3 

A 

1 

2 

3 

4 
1 
2 

3 

4 
1 
2 

3 

4 
1 
2 

3 

4 
1 
2 

3 

4 


N 

0 

0 

0 

0 

1 

1 

1 

1 

2 

2 

2 

2 

3 

3 

3 

3 

4 
4 
4 

4 
3 

5 
5 
3 


10 

to 



REFERENCES 


1. Shanley, F. R. , Simplified Analysis of General Instability of Stiffened 
Shells in Pure Bending, Jour. Aer. Sci. vol. 16, no. 10, pp. 590-592, 
October 1949. 

2. Van der Neut, A. , The General Instability of Stiffened Cylin drical Shells 
under Axial Compression, Rep. S314, National Aeronautical Research 
Institute, Amsterdam, in Reports and Transactions Vol. Xin, 1947, 

pp. S57-S84. 

3. Hedgepeth, J. M. and Hall, D. B. , The Stability of Stiffened Cylinders , 

ER 13731, Space Systems Division, Martin Company, Baltimore, Md. , 
December 1964. 

4. Card, M. F. , Preliminary Results of Compression Tests on Cylinders 
with Eccentric Longitudinal Stiffeners . NASA TM-X-1004, Sept. 1964. 

5. Baruch, M. and Singer, J. , Effect of Eccentricity of Stiffeners on the 
General Instability of Cylindrical Shells under Hydrostatic Pressure, 

Jour. Mech. Eng. Sci., Vol. 5, no. 1, pp. 23-27, March 1963. 

6. Van der Neut, A. , General Instability of Orthogonally S tiffened Cylindri- 
cal Shells, Collected Papers on Instability of Shell Structures - 1962, 

NASA TN-D-1510, pp. 309-321, 1962. 

IT I? 

7. Flugge, W. , Stresses in Shells , Springer-Verlag, Berlin/Gottingen/ 
Heidelberg, 1960. 

8. Almroth, B. O. , Buckling of Orthotropic C ylind ers Under Axial Com pres- 
sion, LMSC-6-90-63-65, Lockheed Missile and Space Company, Sunnyvale, 
California, June, 1963. 

9. Seely, F. B. and Smith, J. O. , Advanced Mechanics of Materials , Second 
Edition, Chapter 6, paragraph 51, John Wiley & Sons, Inc. , New York, 
London, 1952. 

10- IBM 1620 Fortran (wit h format) , File Number 1620-25, Form C26- 
5619-4, IBM Product Publications Department, San Jose, California, 

1963. 


30 


NASA-Langley, 1966 M334 



